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ABSTRACT 

We have studied the opacity of dust grains at submillimeter wavelengths by estimating the optical 
depth from imaging at 160, 250, 350, and 500 /im from the Herschel Gould Belt Survey and comparing 
this to a column density obtained from the 2MASS-derived color excess E(J — K S ). Our main goal was 
to investigate the spatial variations of the opacity due to 'big' grains over a variety of environmental 
conditions and thereby quantify how emission properties of the dust change with column (and volume) 
density. The central and southern areas of the Orion A molecular cloud examined here, with Ah 
ranging from 1.5 xlO 21 cm -2 to 50 x 10 21 cm -2 , are well suited to this approach. We fit the multi- 
frequency Herschel spectral energy distributions (SEDs) of each pixel with a modified blackbody to 
obtain the temperature, T, and optical depth, r 12 oo, at a fiducial frequency of 1200 GHz (250 /im). 
Using a calibration of Ah/ E(J — K s ) for the interstellar medium (ISM) we obtained the opacity (dust 
emission cross-section per H nucleon), cr c (1200), for every pixel. From a value ~ 1 x 10~ 25 cm 2 H _1 at 
the lowest column densities that is typical of the high latitude diffuse ISM, <t c (1200) increases as A^ 28 
over the range studied. This is suggestive of grain evolution. Integrating the SEDs over frequency, we 
also calculated the specific power P (emission power per H) for the big grains. In low column density 
regions where dust clouds are optically thin to the interstellar radiation field (ISRF), P is typically 
3.7 x 10~ 31 W H _1 , again close to that in the high latitude diffuse ISM. However, we find evidence 
for a decrease of P in high column density regions, which would be a natural outcome of attenuation 
of the ISRF that heats the grains, and for localized increases for dust illuminated by nearby stars or 
embedded protostars. 

Subject headings: Dust, extinction - evolution - Infrared: ISM - ISM: structure - Submillimeter: ISM 



1. INTRODUCTION 

Thermal dust emission is optically thin at submillime- 
ter wavelengths. As such, it provides a useful probe of 
the interstellar medium (ISM) and the embedded fila- 
mentary and clumpy structures within it which relate to 
the early stages of star formation. These structures are 
being revealed in exqu isite detail by the Herschel Gould 
Belt Survey (HGBS, lAndre et al.rf2010h . the Herschel 
imaging survey of OB Young Stellar objects (HOBYS, 
iMotte et al.l l2010i ). and the Herschel infrar ed Galactic 
Plane Survey fHi-GAL. iMolinari et al1l2010h . The best 
spatial resolution is obtained by the HGBS given the rel- 
ative proximity of the target molecular clouds. 
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For quantitative measurements of the total column 
density, Ah, and for assessment of the mass and grav- 
itational (in)stability of any structures in the molecu- 
lar cloud, the dust opacity is required. Evidence for 
significant environmental changes in opacity, ranging 
over an order of magnitude , has already been presented; 
see, e.g ., Juvela et al.l (1201 ID: IPlanck Collab oratio n et al.l 



(l2011ailb¥ lMartin et all (120121). Using new HGBS data, 
our goal here is to investigate whether or not there is 
any systematic dependence of opacity on column (and 
volume) density. 

Theoretically, dust grains are expected to evolve in 
a sufficiently dense medium, by aggregation or growth 
of ice mantles. Certainly this is the conclusion of 
a number of numerical simulati ons of dust evolution 
in dense ISM environme nts (e.g.. lOssenkopf fc Henning] 
119941: lOrmel et ail 1201 1[ ). Given the complexity of dust 
in the ISM, however, a completely ab initio prediction 
for a particular environment is challenging, and so it is 
useful to seek observational constraints on the onset and 
magnitude of any environmental changes. 

Volume density, not column density per se, ought to 
be the determinant of dust evolution. For this study we 
have selected the southern part of the Orion A molecu- 
lar cloud imaged by the HGBS (Polychroni et al. 2012 
in preparation) . Its high column density, despite its rela- 
tively high Galactic latitude, b ~ —19°, suggests that the 
emission comes largely from a single cloud, in which case 
the high column density can be reasonably attributed to 
high volume density. Also, at least to a first approxima- 
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tion, the cloud should be bathed in a relatively uniform 
interstellar radiation field (ISRF). For this sort of analy- 
sis, such a field is therefore much more favorable than one 
in the Galactic plane where a wide range of conditions 
would be superimposed along the line of sight. 

The larger ('big') dust grains, which account for most 
of the dust mass, are i n thermal equilibrium with the 
ambient radiation field ( Compie gne et ai1l2011l) . The in- 
tensity of the thermal dust emission, when optically thin, 
is given by 

I v = t v B v (T) = a e {v)NuB v {T), (1) 

where t„ is the dust optical depth of the column of mate- 
rial, B V {T) is the Planck function for dust temperature 
T, Nh is the total hydrogen column density (H in any 
form), and u c {v) ~ t^/A^h is the opacity of the interstel- 
lar material (the emission cross-section per H nucleon) F 1 ! 

Clearly, quantifying the optical depth and the opac- 
ity requires knowledge of T. Determination of T is 
made possible by the multi-frequency cove rage now avail- 
able w ith the Herschel Space Observatory (jPilbratt et al.1 
2010), through fitting the spectral energy distribution 
(SED) of I v . The other requirement is an independent 
measure of Nr. Here we have used the near- infrared 
color excess E{J — K$) as a proxy (see Appendix IX)) , de- 
rived from the Two Micron All Sky Survey (2MASS3) 
data. 

Another quantity that we examine (Section [5j is the 
'specific power,' i.e., the total power per H emitted by 
dust grains in thermal equilibrium: 

P = J ^o c (v)B v {T)dv. (2) 

Since in equilibrium P is equal to the total energy ab- 
sorbed per H, P is affected by the intensity of the ISRF, 
which can be higher than average near a strong source of 
radiation, or lower than average because of attenuation in 
a region of high column density (without internal sources 
of illumination). Note how, for a given absorbed P, the 
resulting equilibrium temperature is inversely related to 
how efficiently grains can emit. 

This paper is organized as follows. In Section [2] we 
briefly describe Herschel ima ging of the Orion A m olec- 
ular clo ud using the PAC S (|Poglitsch et al.l 120101 ) and 
SPIRE (|Griffin et alj|2010h cameras. Maps of the SED 
fitting parameters r and T are derived in Section |3J We 
discuss the SED fitting, validate the fits through predic- 
tion of the 100 /im emission for comparison with observa- 
tions by IRAS (Appendix IB1 . and assess the uncertainty 
in /„, and estimate the uncertainties of the derived pa- 
rameters through Monte-Carlo simulation (Appendix [C]). 
In Section 21 we compare r with A^h to find an estimate 
of the opacity (Equation (JTJ). It appears that the opacity 

11 Note the correspondence a c (v) = nmnrK v where re„ is the 
mass absorption (or emission) cross-section per gram of dust, r 
is dust-to-gas mass ratio, and p, is the mean weight per H (1.4). 
Only the product tk v is needed; this quantity is often also called 
the opacity, now in the alternate units cm 2 g — 1 . 

12 The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 



grows systematically with A%, evidence for grain evolu- 
tion. We examine the dependence of quantities related to 
P in Section [5] specifically the anti-correlation between 
T and r and the spatial and Nh dependence of P. This 
analysis provides insight into the various interrelation- 
ships between T, a c , and P discussed in Section [B] We 
conclude with a summary in Section [7] 

2. HERSCHEL OBSERVATIONS 

As part of the HGBS, three fields in the Orion A molec- 
ular cloud were mapped at a scanning speed of 60" s _1 
in parallel mode, acqu iring data simultane ously in five 
bands using th e PACS (|Poglitsch et al.|[2010h and SPIRE 
( Grif fin et al.l IMch arrays. Images were produce d by 
the ROMAGAL map-maker (jTraficante et al.H201lD and 
first results on the filamentary and core substructures 
in the two fields studied here, called Orion A-Cl and 
Orion A-Sl, are presented by Polychroni et al. (2012). 
The images have angular resolutions of 9'.'6, 13'.'5, 18'.'0, 
24'.'0, and 37"0 at 70, 160, 250, 350, and 500 /zm, re- 
spectively. In the following analysis we did not use the 
70 /zm image. Zero-point offsets added to these images, 
obtained by correlat ing with Planck and IRAS images 
(|Bernard et al.ll2010t ), were 8, 20, 10, and 4 MJy sr^ 1 at 
160, 250, 350, and 500 /zm, respectively for the Orion A- 
Cl images and similarly 3, 14, 18, and 8 MJy sr _1 for 
the Orion A-Sl images. Uncertainties in the zero-point 
offsets have the most effect on the SEDs of pixels with 
lower brightnesses, and the offsets were refined as dis- 
cussed in Appendix [Cj Prior to fitting SEDs for each 
pixel, we convolved and then regridded the individual 
images to correspond to the lowest resolution (37'.'0) on 
a common grid with 11'.' 5 pixels. Power spectra of these 
images follow a 'cirrus-like' power-law relation that de- 
cays with spatial frequency until at the highest frequen- 
cies it merg es with a flat level corresponding to the noise 
in the map (|Rov et al.ll2010HMartin et al.l l2010), assessed 
at 0.92, 0.58, 0.31, and 0.38 MJy sr" 1 for 160, 250, 350, 
and 500 /zm, respectively. 

3. MAPS OF r AND T 

We parametrized the spectral dependence of the opac- 
ity (or t) as a c (v) / <j c (vq) = (v/vq)?, with a fidu- 
cial frequency Vq = 1200 GHz (250 /zm). Although 
HGBS consortium assumed an emissivity index, /?, 
of 2.0, however, in this paper we adopted a fixed 
= 1.8 to facilitate comparison with earlier analyse s 
(jPlanck Collaboration et al.M2011at IMartin et all 1201 2j) . 
Treating /3 as a fitting parameter does not change our 
results (Section |4.1[) . The SEDs are thus fit with two pa- 
rameters, T1200 and T. In our SED fits, we have not used 
70 /zm PACS data because of probable contamination due 
to non-equilibrium emission by smaller dust grains (Very 
Small Grains, or VSGs) which broadens the apparent 
SED toward wavelengths shortward of the spectral peak. 
The SED of cold dust emission for j3 = 1.8 and a temper- 
ature of 15 K, typical of this cloud, peaks at 200 /im, and 
so the remaining four submillimeter Herschel passbands 
are still sufficient to undertake this study. The data fit 
this simple model of a modified blackbody function well. 
A representative SED and its fit are shown in Figure [TJ 

The IRAS 100 /zm brightness was used in fitting SEDs 
in previous work (|Pfanck Collaboration et all 12011a; 
IMartin et al.l l2012h . but was not included here because 
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Figure 1. Representative SED in the Orion A map. The derived 
temperature is 17.9 K and the optical depth at 250 (im is 6.3xl0 — 4 . 
Only the four lowest frequency Herschel passbands were used in fit- 
ting the modified blackbody SED. This model predicts the 100 fim 
IRAS brightness quite well (open circle, see Appendix [EJ but un- 
derpredicts the 70 £im emission (triangle) where non-equilibrium 
emission from VSGs becomes important. 

of its relatively low angular resolution (~ 4.'3). As dis- 
cussed in Appendix [Bj we subsequently checked that the 
100 /im emission is consistent with the SED fit. On the 
other hand, the 70 /im emission is greater than that pre- 
dicted by the equilibrium emission SED, confirming the 
additional contribution from VSGs (see example in Fig- 
ureH}. 

Errors on the parameters from the SED fit were esti- 
mated using Monte Carlo simulations after assessing the 
various sources of error in the intensities of the Herschel 
images (see Appendix [C]). A typical SED fit has a % 2 
near two (as expected for four intensities, two parame- 
ters, and thus two degrees of freedom) and the map of 
X 2 is featureless. 

Figure [2] shows results from the SED fits as maps of 
the fitting parameters for the Orion A-Cl field. The 
left panel shows a map of T1200; the variation of Ah is 
shown by the top color bar, adopting the mean relation 
found below (Figure [3]). The corresponding map of T is 
in the right panel; a complementary inverted logarith- 
mic scale has been used to bring out both the general 
anti-correlation (see also Figure [5] in Section [Oil and ex- 
ceptions. One common feature in these two maps of the 
fitting parameters is the filamentary structure, which is 
of high column density and generally low temperature, 
as might be expected for an IRSF that is attenuated in 
those regions. This behavior is discussed fu rther in Sec 
tionEl Filaments are comm o n in th e HGBS (lAndre et al 



201C ; iMen'shchikoy et all 120101: lArzoumanian et al 
20rl iPeretto et all 120121). HOBYS dMotte et al.ll2010t 



Hill et al.ll201lHHennemann et alj|2012t ISchneider et alj 
2012J), and Hi-GAL ( jMolinari et al.l l2010[) images. 



4. SUBMILLIMETER OPTICAL DEPTH, JV H , AND 
OPACITY 

4.1. r and A% 

Before comparison with our E( J — K s ) proxy for Ah, 
all parameters resulting from the SED fits on the Her- 
schel maps were sampled to match the map of E(J — K s ) 
(Appendix IA.3|) . Figure [3] shows a plot of the sampled 
r i200 agai nst A h from the conversion from E(J — K s ) in 
Equation (|A1|) . 

The two independent measures of column density are 



well correlated. The data have some dispersion, with pos- 
sible reasons including uncertainty in the measurement 
of color excess, cosmic scatter about Equation (|A1I) con- 
verting E( J — K s ) to Ah, and errors in the determination 
of T1200 from the SED. 

Obscrvationally, the properties of dust grains have 
been best constrained for the high latitude dif- 
fuse ISM. For example, the equilibrium dust tem- 
perature is typically near 18 K and a e (1200) 
is 1.0 x !Q- 25 cm- 2 H" 1 (TBoulanger et all 119961: 
iPlanck Collaboration et all 1201 laf) . The expected lin- 
ear relation for this opacity is plotted with the dotted- 
dashed line in Figure [3J This relation is close to that 
of the present data at low column density. However, 
there appears to be a significant non-linearity in the data. 
The formal fit, excluding Ah data with S/N < 2, is 
T1200 7V H ' 28±0 01 ; the systematic error on the index is 

0. 03. Such non-linearity appears to be evidence for grain 
evolution. 

Recall that T1200 was derived under the assumption of a 
constant (3. We also derived r 12 oo treating (3 as a (third) 
free parameter, in case /? evolves too. However, we found 
that j3 does not change systematically with column den- 
sity, and so the same non-linear trend is observed in Fig- 
ure[3]- the evidence for grain evolution remains. There is 
increased scatter about the mean correlation and the in- 
tercept is 10% higher because the average f3 is 1.96 (with 
higher (3, SED fits tend to have lower T and hence higher 
r). 

Although the conversion factor from E{J — K s ) to 
Ah might not be constant for the higher Ah in Fig- 
ure |31 there is already some divergent non-linear behavior 
within the range A h< 5 x 10 21 cm -2 that is calibrated 
(|Martin et al.l l2012f) . Nevertheless, we have to be aware 
that apparent non-linearity might simply be a reflex of 
an er rone ous conversion of E{ J — K s ) to N-r (see Equa- 
tion (|A1[1 ). To remove the evidence for grain evolution, 

1. e., to maintain a constant submillimeter opacity, would 
require Ah oc E(J — A s ) 128 - which would itself be ev- 
idence for grain evolution, albeit as it affects the near- 
infrared extinction. This dependence, is however, oppo- 
site to what has been predicted for the initial changes in 
extinction curves resulting from gr ain evolution by i ce- 
mantle formation and aggregation (jOrmel et al.ll201l"l) . 

Furthermore, empirically, the ratio of E(J "— H) / ' E(H - 



A s ) is constant wit h colu mn density (Appendix IA.1|) 



I Chapman fe Mundv (2009) find the same constancy in 
the shape of the near-infrared extinction curve in probes 
of cores to slightly higher peak optical depth than 
reached here, but interestingly they find an increase in 
the relative amount of mid-infrared extinction at high 
optical depth. A constant shape of the near-infrared ex- 
tinction curve, despite changes in the visible (e.g., as en- 
coded by Rv), would be consistent with aggregation of 
the smallest interstellar grains with the largest, but not 
with aggregation among the largest particles, because 
in the latter case the scattering would be greatly en- 
hanced contribut ing to a steeper wavelength dependence 
(|Kim et al.l 119941 ) . The enhanced mid- infrared extinc- 
tion seen at high column den sities might be explained b y 
the addition of ice mantles (Chapm an fc Mundvl 120091) . 
Possibly the same affects the submillimeter. Grain evo- 
lution, when and where it occurs, might be complex and 
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Figure 2. Left: Optical depth map for Orion A-Cl as obtained from fitting SEDs to the 160-/im PACS and three SPIRE intensities, from 
images brought to a common resolution of 37'.'0. The dust emissivity index fj was fixed at 1.8. Right: Corresponding temperature map, 
with a complementary inverted color scale. The logarithmic range spans 12.3 to 20.0 K. Note by comparison with the optical depth how 
the dust in regions of high column density tends to be cooler (< 1.14 dex or 13.8 K) than in diffuse low column density regions (> 1.23 dex 
or 17 K). Within the cool filaments are a few compact regions where the dust is heated by embedded protostars. The few more dif fuse 
regions with low column density but even higher dust temperature have an enhanced local radiation field; sec also Figure [5] in Section 15.21 
below. 
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different than this and so it is hard to be definitive. Until 
there is evidence to the contrary, we think it is reason- 
able to adopt E(J — K s ) as a measure of Nn over the 
range of column densities found in this field, and thus 
we conclude that it is the submillimeter opacity that has 
changed. But it remains important to understand how 
the submillimeter opacity can change without a notice- 
able change in the near infrared. 

Dense filaments appear to have a charac t eristic 
linear width ~ 0.1 p c (|Arzoumanian etai] l20ll 
iFischera fc Martid 12012$ ) which is 50" at the ~ 400 pc 
distance of Orion. This scale is not much larger than 
the 37" resolution of our Herschel maps of T1200 but 
is significantly smaller tha n the ~ 3'5 resolution of the 
E(J ~ K s ) map ( Appendix IA. 21) required for our opacity 
study. Therefore, we do not have enough spatial reso- 
lution to investigate in detail the opacity of dust inside 
individual filaments (or prestellar cores) where the col- 
umn density and volume density would be even higher. 
We are not able to see if the trend indicating grain evolu- 
tion continues to even higher densities. But see the brief 
discussion ending the next subsection. 

The highest column densities (Ay > 8) are the most 
likely to have been influenced by self gravity. In prin- 
ciple, if extremely compact structures developed there 
and had a low covering factor they might be missed 
in the 'AvMAP' assessment of E(J — K s ) because of 
the rela tivel y sparse sampling by the 2MASS stars (Ap- 
pendix IA.1I) . and yet they would still contribute to 
the Herschel assessment of T1200 because their emission 
would be included in the finite beam. If that were the 
case and if such structures had a significant column den- 
sity, then the Ah derived from E(J — K s ) would be un- 
derestimated relative to the sampled T1200, and so the 
observed trend in Figure [3] would turn up above this high 
column density. If anything, the data there fall slightly 
below the non-linear dependence shown. This mapped 
region of the molecular cloud does have some com- 
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Figure 3. T1200 from SED fits versus column density Nn ob- 
tained from near-infrared color excess E(J — K B ). Given the large 
number of pixels in the maps, we have used a logarithmic grey- 
scale to represent different values in a two-dimensional histogram 
of the number of pixels per bin. Individual data points are plot- 
ted (dots) where the density of points is low, less than or equal to 
three. Data where Nn has S/N < 2 are plotted as red pluses and 
are not included in the histogram or fit; these are also identified in 
subsequent figures. The crosses show data binned logarithmically 
along Nft. The dotted-dashed line shows a linear relationship be- 
tween T1200 and the underlying TVh based on a fixed dust opacity, 
here <r c (1200) = 1 X 10~ 25 enfi IT" 1 typical of high latitude dust. 
Compared to this relation there appears to be a non-linear increase 
in T1200 (and thus <r e (1200)) with JVjj- The solid line is the power- 
law correlation obtained from fitting r and 7Vh, excluding the low 
S/N data. The fit yields a power law index 1.28 ± 0.01 (slope of 
the line) and amplitude (1.01 ±0.02) X 10 — 4 for the adopted units 



(intercept at JVh = 10 cm ). 
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Figure 4. Opacity <r c (1200) versus Nn. The solid line shows a 
power-law variation <r c (1200) oc Nh 0,2s taken from the solid line in 
Figure [3] The dashed line depicts the effect of color correction due 
to the finite bandwidth of the 2MASS filters on the measurement 
of E(J — K B ) (column density). The horizontal dotted-dashed line 
shows the average <r c (1200) for the diffuse ISM. The dotted line 
is for <7 C (1200) = 3.3 X 10" 25 cm 2 H" 1 , the standard adopted 
by the HGBS and HOBYS. Pixels representing two-dimensional 
histogram is same as described in Figure [3] Individual data points 
are plotted where the density of points is less than or equal to 
three. 

pact (self-gravitating and probably quasi-equilibrium) 
high column density structures, filaments and even some 
'cores,' with enough mass to be detectable in the sub- 
millimeter, but their physical size is large enough that at 
the distance to the Orion A cloud these structures are 
measured to be resolved, and so are not as in our hy- 
pothetical scenario physically vanishingly small and not 
sampled by the 2MASS stars. Another constraint on 
any hypothesis aiming to account for, or explain away, 
this non-linearity is that the trend in Figure [3] begins at 
quite low column densities. Furthermore, there is ample 
evidence for a higher opacity associated with high col- 
umn density mole cular regions (see, e.g., the summary 
in Appendix A in iMartin et al.l 120121 and the following 
subsection). 

4.2. Opacity <j c (1200) and N K 

The nature of the systematic change in a c (1200) is re- 
vealed more explicitly in Figure |4j which is just an al- 
ternative representation of Figure |31 from the sampled 
^1200 map simply divided by the JVh map, plotted on 
an expanded vertical scale. The solid and dotted-dashed 
lines are the same as defined for Figure [3] Despite the 
dispersion about the solid line, the overall trend clearly 
indicates an increase of <7 e (1200) with column density 
(cr c (1200) oc jV°- 28±0 - 01±0 - 03 ), by at least a factor of two 
over this Ah range. 

A particular additional perspective provided by the 
present study of this region in Orion A is that there 
ought to be a rough correspondence between volume 
density and column density. This link could enable a 
more direct connection of grain evolution to environment. 
Compared to the density of the diffuse ISM (nn ~ 1 
10 cm -3 ), the density typical of a star formation region 
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where clumps are beginning to show signs of gravitational 
contraction is several orders of magnitude higher (10 3 - 
10 5 cm" 3 ). In the present analysis, we are probing envi- 
ronments, especially in filamentary structures, with «h 
up to ~ 10 4 cm" 3 . 

The measurement of er c (1200) is sensitive to the esti- 
mate of temperature through the Planck function; see 
Equation (fTJ) where we have assumed a constant T along 
the line of sight. Obviously, either for passively-heated 
molecular clouds with significant optical attenuation or 
for those with an internal energy source, there is a gra- 
dient in temperature which we have ignored. Particu- 
larly relevant here is the former possibility which has 
been studied extensi vely using ra d iative transfer simula- 
tions. For example, lYsard et al.l (|2012[ ) found that the 
temperature obtained from fitting a single-temperature 
SED is overestimated compared to the true dust tem- 
perature for central lines of sight through filaments or 
spherical clouds. The difference increases non-linearly 
with central column density and is noticeable by Ah ~ 
14 x 10 21 cm" 2 , depending on the grain model. In the 
context of our analysis, an 'overestimate' of the temper- 
ature would have suppressed the optical depth derived 
here; hence, this effect would enhance rather than erase 
the trend of increasing ct o (1200) with Ah- 

The dashed line in Figure |4] shows the correction of 
the correlation line when the effect of the finite filter 
bandwidths on the color excess is taken into account (see 
Appendix [A3])- Underestimation of E(J — K S ) makes the 
rate of increase in a e (1200) only marginally smaller. 

For the pixels with low column density, the de- 
rived values of a c (1200) are close to the value 1 x 
10" 25 c m 2 H" 1 for high lati tude diffuse interstel- 
lar dust (jBoulanger et al.l |l996) and typically within 
the range of (0.6 — 1.6) x 1 0~ 25 cm 2 H" 1 found 
by IPlanck Collaboration et al.l (|2011aD . We find that 
er c (1200) increases systematically with Ah by at least 
a factor two. Such high opacities are not unprece- 
dented, one example being 3.8 x 10 -25 cm 2 H" 1 for 
a dense molecular region in the Vela molecular ridge 
where the column density r anged between 10 to 40 x 10 
cm" 2 (jMartin et al.l [20121) . For a molecular region in 
Taurus of intermediate column den sity (~ 10 22 cm" 2 ), 
IPlanck Co llaborati on et all (|2011b|) obtained an opacity 
of (2.3 ± 0.3) x 10" 25 cm 2 H -1 . By contrast, for the 
neighboring atomic phase with Ah < 3 x 10 21 cm -2 , they 
obtained (1.14±0.2) x 10~ 25 cm 2 H" 1 . Based on a com- 
plementary new method, physical modeling of the br ight- 
ness profiles of filaments. iFischera fc Martini (|2012af) find 
a range (1.3 — 2.8) x 10" 25 cm 2 H _1 (the range depending 
on the adopted distance). These are consistent with the 
trend seen in Figure (?] 

The standar d opacity adopte d by the HGBS a nd 
HOBYS (e.g., lAndre et all l20loT: IMotte et al.l l20Toh is 
cr e (1200) = 3.3 x 10" 25 cm 2 H" 1 (see Figure H). This 
is take n from theoretical calculations by Prcibisc h et al.l 
(1 19931) for evolved dust in protostellar cores (see also 
IQssenkopf fc Henning| [T994). Although we do not probe 
to such high densities, this value seems a reasonable ex- 
trapolation of the trend in Figures [3] and [U a corollary is 
that a single value of the opacity cannot be used for all 
environments across a region. 

As emphasized by [Shirley et al.l (|2011f ). the largest un- 



6 



Roy et al. 



certainty in determining the mass of a core is in the 
adopted value of the opacity. Their modeling of the 
extended submillimeter emission in the cold envelope 
of B335, combined with deep imaging at H and K 
{E(H — K) up to 3, about three times larger than in 
our map), allowed them to constrain the ratio of the 
submillimeter opacity, at 850 /im, to the opacity (scat- 
tering plus absorption) at 2.2 /im. They found a range 
(3.2 — 4.8) x 10~ 4 . Extrapolating this to the ratio in- 
volving the opacity at 250 fiin using their estimate of /? 
(2.18 to 2.58) gives a range (4.6 - 11.3) x 10~ 3 . This 
can be compared to our estimate of the same ratio, 
(2.3-4.9) xl0~ 3 for cr c (1200) = (l-2)xl0~ 25 cm 2 H" 1 
(note that this ratio is independent of the calibration of 
either opacity to Nr). Thus regions characterized by 
higher density, or higher column density here, both ap- 
pear to have an enhanced submillimeter opacity relative 
to that in the K band (and one cannot rule out that both 
might have increased). Of course, as we have noted in 
the Introduction, what is really needed to find the mass is 
the product rn v , and so without some direct calibration 
as attempted here, the Achilles heel is to first derive the 
submillimeter dust opacity itself from the above opacity 
ratio and then adopt an appropriate dust to gas ratio r. 

5. POWER 

Big dust grains bathed in the ISRF absorb energy at 
higher frequencies, and in thermal equilibrium they re- 
emit the same amount of energy at much lower frequen- 
cies. Thus, the observed P is what the grains were able 
to absorb in their environment, per H. For the dust emis- 
sivity relation (/3) adopted, P can be expressed as 

P/P = (o- a (1200)/o-o) (T/T ) 5 - 8 , (3) 

normalized in terms of the average paramc- 
ters in the diffuse atom i c ISM at high latitude 
(jPlanck Collaboration et~aH l2011al ): T = 17.9 K, 
a = 1.0 x 10" 25 cm 2 H" 1 and P = 3.8 x 10" 31 W H _1 . 

5.1. T and r 

The value of the equilibrium grain temperature is not a 
primary dust parameter. Rather, as is clear from Equa- 
tion (|3]), T is a parameter that responds to environmen- 
tal changes in the ISRF affecting the absorbed P and to 
changes in the grain properties like ct c (1200) and possibly 
the absorption which also affects the absorbed P. 

The anti-correlation of T and T1200 seen in the maps of 
Figure [5] and similarly for the Orion A-Sl field is high- 
lighted in Figure [5j Because we are primarily interested 
in regions where there is no internal source of energy, we 
have not plotted data for the pixels corresponding to a 
few identified embedded protostars (high T and interme- 
diate Ti2oo)- The dashed curves are not a fit but rather 
trace the general trend schematically. We have trans- 
ferred these loci to Figures [7] - [9] below to gain insight 
into relationships between other variables. For simplicity, 
we have used two power laws, T/Tb = (T(1200)/rb)~ a , 
where a break occurs at Tb = 14.7 K and Tb = 4.9 x 10~ 3 , 
and where a — 0.06 and 0.1 below and above the break, 
respectively. 

5.2. Spatial Dependence of P 
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Figure 5. T versus T1200 from the SED fits, with the logarithmic 
grey-scale representing the values in a two-dimensional histogram. 
Representative error ellipses for a single pixel are from Monte Carlo 
simulations of the SED fitting. The dashed curves are power laws 
schematically following the anti-correlation. Individual data points 
are plotted where the density of points is less than or equal to three. 
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Figure 6. Map of the specific power P in the Orion A-Cl 
field. The overplotted white contour, corresponding to TVjj = 
20 X 10 21 cm~ 2 (Ay ~ 10.6 mag), reveals how P decreases sig- 
nificantly in high column density regions, especially here along fil- 
aments. 



We produced a map of P shown in Figure [6] making 
use of the T and r maps (as in Figure^ and the average 
dependence of er c (1200) on T1200 from Figures [3] and 2J 
More directly from the data, a map of the radiated energy 
is obtained from the integration of the respective SEDs 
over all frequencies. As in Section 14.21 we divide the 
resulting sampled map by the Nn map to obtain a map of 
the specific power P (Equation ([2])). This is very similar 
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to Figure E0 

For the pixels corresponding to low Ah, P is on average 
about 3.7xl0~ 31 W H _1 , close to the value 3.8 x KT 31 
W H" 1 for the typical speci fic power in the diffuse ISM 
at hig h Galactic latitude ([Planck Collaboration et all 
I2011al> . In a few diffuse regions in our maps, for ex- 
ample the northern part of Orion A-Cl, P gets as high 
as 10 x 1CP 31 W H _1 , clearly showing an enhancement 
of the local radiation field by a factor of ~ 2.5 times rel- 
ative to the average ISRF. Others examples are the B77 
Bright Nebula and a halo around V1792 Ori, both in the 
Orion A-Sl field and so not shown here. These regions 
also have strong emission at 70 /urn and 160 /im in the 
respective Herschel images. 

As illustrated in Figure [6] by the overlay of the contour 
Ah = 20 x 10 21 cm -2 , regions with the highest Ah tend 
to have lower P. The most straightforward qualitative 
interpretation is in terms of attenuation of the ISRF. The 
high frequency components of the ISRF that are primar- 
ily responsible for the heating of the dust grains cannot 
penetrate easily into the interior of the denser molec- 
ular structures. Most of the power is absorbed within 
the outer layers of the cloud and so the dust in the core 
generally has a lower equilibrium temperature. 

It should be noted that only a fraction of the total Ah 
through the center of the filaments is associated directly 
with the structure, the rest arising from the embedding 
and foreground and background material. Therefore, the 
detailed relationship between specific power and column 
density requires knowledge of density profiles and cloud 
geometry, the strength of the ambient ISRF, and detailed 
radiative transfer modeling taking into account appro- 
priate dust populations and absorption an d scattering 
properties (|Fischerall20TH lYsard et al.ll2012t ). 

In contrast, and not surprisingly, close inspection of 
Figure [6] reveals a few lines of sight where there is an 
internal source of energy (protostar(s)) within a dense 
compact structure. With this extra energy input, the 
resulting P is higher locally. These are readily seen in 
Figure[2] (right) as warmer than the immediate surround- 
ings. An example is at RA, Dec 05:39:56, —7:30:27, in 
a complex environment with the brightest hot source 
named LDN 1641 S3 IRS/ FIRSSE 101. 

5.3. P and A H 

Figure [7] shows the relationship of sampled P to Ah- 
There is a trend of decreasing P at higher column den- 
sity that sets in well before the value 20 x 10 21 cm" 2 
emphasized in Figure [5] Above we quantified trends 
relating c c (1200) and T1200 to Ah, and T to ti 2 oo- If 
we substitute these trends in Equation ([3]), we find 
P/P h = (N H /Nn h )~ s , where P h = 2.98 x 10" 31 W IT 1 
and A Hb = 20.3 x 10 21 cm" 2 , and where S = 0.165 and 
0.46 below and above the break, respectively. These loci 
are plotted in Figure [71 

6. INTERRELATIONSHIPS AMONG T, P, AND o- e (1200) 

Among T, P, and ct c (1200), a scatter plot of any two 
provides insight into their interrelationship. On the re- 

13 The displayed map of P is shown at higher spatial resolution 
than a map of <r e (1200) would have to highlight the effects of at- 
tenuation of high frequency radiation in the filamentary structures. 
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Figure 7. Specific power P versus column density ATh- P fluctu- 
ates about the typical high latitude diffuse ISM value of 3.8 X 10~ 31 
W H (dotted-dashed line). For high column density, P decreases 
because of attenuation of the ISRF. Dashed lines show the loci 
transferred combining trends from Figures [3] - E] Pixels represent- 
ing two-dimensional histogram is same as described in Figure [3] 
Individual data points are plotted where the density of points is 
less than or equal to three. 
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Figure 8. Dust temperature T versus dust opacity tT e (1200). 
Solid lines are the loci of constant specific power, P. The middle 
line is for 3.8 xl0~ 31 W H" 1 , the typical P in the high latitude 
mean ISRF (see Section [5J. The lower and upper lines are for an 
ISRF different by factors of 0.5 and 1.5, respectively. Dashed lines 
show the loci transferred combining trends from Figures[3]-[5] Pix- 
els representing two-dimensional histogram is same as described in 
Figure [3] Individual data points are plotted where the density of 
points is less than or equal to three. 

suiting figures, we also plot loci defined by fixing the 
third quantity. 

6.1. T - a c (1200) Relation 

As discussed above, the dust temperature is deter- 
mined by the energy balance between the power absorbed 
from the ISRF and emission. For a given intensity of 
the ISRF and absorption coefficient, and thus absorbed 
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Figure 9. Specific power P versus dust opacity <r c (1200). The 
solid lines (bottom to top) are the loci of constant temperatures 
equal to 13, 14, 15, and 17.9 K. Dashed lines show the loci trans- 
ferred combining trends from Figures \3\ - \5\ Pixels representing 
two-dimensional histogram is same as described in Figure [3] Indi- 
vidual data points are plotted where the density of points is less 
than or equal to three. 

power, grains with a higher cr c (1200) radiate more ef- 
ficiently and so attain a lower equilibrium temperature 
(Equation [3]) . Therefore, the observed equilibrium dust 
temperature and the intrinsic dust property ov-f^OPf 1 "^ 
are expected to be inversely related. In Figure[8l we show 
the anti-correlation between sampled T and <t c (1200) for 
the Orion A region. The solid line shows a locus of con- 
stant specific power P — 3.8 x 10~ 31 W H _1 whereas the 
loci below and above are 0.5 and 1.5 times this value, 
respectively. There is a range of P and indeed, unlike in 
the high latitude diffuse ISM, we do not expect to have 
constant P in the regions of high column and volume den- 
sity (see Figure [7]), due to attenuation by dust itself. To 
capture this effect, we have transferred the loci from the 
above figures. We find T/T h = {a e {l2QQ)/a b )-\ where 
cr b = 2.41 x 1(T 25 cm 2 H" 1 and where e = 0.27 and 0.46 
below and above the break, respectively. These loci are 
plotted in Figure [5] 

6.2. P - a c (1200) Relation 

Figure [5] shows sampled P versus <r e (1200), and is suffi- 
cient to complete our examination of the three-parameter 
relations of Equation §5§ (P versus T being redundant). 
For a given dust equilibrium temperature, the power ra- 
diated in the submillimeter is proportional to the emis- 
sion cross-section. The solid lines are the loci of constant 
temperature with increasing slope corresponding to 13, 
14, 15, and 17.9 K, respectively. For dust grains exposed 
to same ISRF (constant P), the less efficient emitters 
(lower <7 o (1200)) will have higher temperature. 

To capture the effect of attenuation, we have trans- 
ferred the loci from the previous figures. We find P/Pb = 
((Te(1200)/crb) _t , where i — 0.59 and 1.65 below and 
above the break, respectively. These loci are plotted in 
Figure [HI 

14 <T e (1200) also involves the dust to gas ratio. 



7. CONCLUSION 

We have studied the properties of dust grains using 
multi- wavelength images of thermal dust emission in the 
Orion A molecular cloud at 160, 250, 350, and 500 /xm, 
acquired by the PACS and SPIRE cameras on Herschel 
as part of the HGBS. We fit a modified blackbody model 
to the spectral dependence of the surface brightness of 
each pixel. Although dust along the line of sight is prob- 
ably not all at the same temperature, assuming a single 
T is a reasonable model providing good quality SED fits 
to the Herschel data. Thus, we obtained the spatial dis- 
tribution of T and T1200 across the map. 

The derived optical depth map is well correlated with 
the iVn map derived from E(J — K s ) and we combined 
them to explore how the dust opacity a e (1200) varies 
with column density. A submillimeter opacity er c (1200) 
near 1 x 10~ 25 cm 2 H _1 is typical for diffuse dust along 
high latitude lines of sight of low column density (A\r~ 
1 mag) and close to that found at low column density 
in these Orion A fields. On lines of sight intercept- 
ing high column density filamentary structures, how- 
ever, we found evidence for an increase of cr e (1200) by 
a factor two or more, approaching that adopted by the 
HGBS and HOBYS. Overall, there is a systematic trend 
er o (1200) cx N^ 2S . This dependence is strong evidence 
for grain evolution with environment. There is not a sin- 
gle opacity that can be applied to a whole mapped region 
with environments ranging from diffuse to higher column 
density structures. This has quantitative implications for 
the interpretation of the probability distribution function 
of the column density and possibly for the shape of the 
core mass function. 

For the low column density lines of sight, the average 
value for the specific power P is 3.7 x 10 -31 W H _1 , 
equivalent to 1.2 Lq/Mq. The power emitted is equal to 
the power absorbed and so depends on the local ambient 
ISRF. Thus, the decrease of P that we see in the dense 
filamentary structures (P ~ 1.5 x 10~ 31 W H _1 ) can 
be attributed to attenuation of the ISRF. There are also 
local enhancements of the ISRF. Overall, P ranges over 
a factor of 10, roughly centered on the diffuse ISM value. 

The emission opacity cr c (1200) of big grains together 
with the relative strength of the ISRF determines the 
equilibrium temperature. In the diffuse ISM where dust 
grains are exposed fully to the ISRF, the temperature 
is typically 18 K. Inside dense filamentary structures in 
the observed region of the Orion A molecular cloud, the 
opacity is larger and the ISRF is attenuated, both leading 
to dust temperatures lower than 14 K. 

In this high latitude region in the Orion A molecu- 
lar cloud, high column density arguably corresponds to 
high volume density. The observed change in the opti- 
cal properties of dust grains in dense cold regions can 
be attributed to grain growth due to aggregation and/or 
accretion of ice mantles. This process might reasonably 
be expected to correlate with a decrease in the relative 
abundance of VSGs. It is difficult, however, to probe the 
presence of the latter because the radiation that would 
normally excite their non-equilibrium emission at shorter 
wavelengths (< 100 /zm) is sharply attenuated in these 
very regions. 

We thank the referee for helpful comments re clarifica- 
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APPENDIX 

A. 2 MASS COLOR EXCESS AND COLUMN DENSITY 

Column density can be estimated by measuring extinction or color excess. Very productive use has been made 
of the near-infrared 2MASS data in the J, H, and K s passba nds, post-processed by a variety of techniques (e.g., 
iRowles Sz Froebrichl 120091 : iDobashil 120 lit lLombardi et al.l 120111 ) . To stay close to the data and avoid unnecessary 
additional assumptions (see also Appendix IA.4I) , here we quantify the extinction in terms of the near-infrared color 
excess E(J — K s ), rather than the total visual extinction Ay. For comparison to values cited in other literature, the 
adopted conversion has been Ay = 5.89 E( J — K s ). 

A.l. Sampling the Color Excess 

The extinction map was derived using the 'AvMAP' procedure (Schneid er et al.ll2011l ) as a weighted average of 
measures of E(J — H) and E(H — K s ), and subsequently expressed here as E(J — K s ). The map is created from a 
weighted mean of the color excesses of individual stars. The weighting function is a Gaussian of FWHM — 4' centered 
on each 2' pixel and stars are considered out to a radius of 3 a (1.3 FWHM). To monitor the quality, maps of the 
weighted number of stars used and the weighted standard deviation of the mean (i.e., the uncertainty of E(J — K s )) 
are also computed. 

Given the high column densities reached in this field, the issue arises whether the color excess maps might somehow 
be 'saturated,' i.e., missing the high est values. This is s ue do es no t seem to be ad dressed in the literature. We have 
made scatter plots of the maps of IRowles fc Froebrichl (|2009f ) and[Dobashi (1201 If ) compared to our E(J — K s ) map. 
These plots show good correlations with the expected slope up to Ay ~ 8, after which these other products cease to 
continue to rise as our E(J — K s ) grows further. 

Maps of E(J — H) and E(H — K s ) can be der ived separately. Th ese co rrelate very well with a slope of 1.73 ± 0.01, 
close to the 'normal' value 1.7 (see discussion in lMartin et al.l l2012). Like lLombardi et al.l (|201lD have also shown for 
the Orion A cloud, this normal color-color relation extends beyond E{J — H) = 1.7 mag and E(H — K s ) = 1 mag, 
though the correlation becomes noisier. Thus, there seems to be no evidence for either the effects of grain evolution 
or of saturation of one color with respect to the other up to E(J — K s ) = 2.7 mag which corresponds to Ay ~ 16 mag. 
Nevertheless, we need to acknowledge the issue and ensure that the colu mn d ensity from Herschel r is sampled on a 
closely comparable basis. This is accomplished as described in Appendix I A. 31 

A. 2. Gauging the Resolution 

In our Herschel maps of T1200 at 37" resolution there is structure on smaller scales than the 4' weighting function 
of 'AvMAP.' Therefore, it is important that the map of T1200 be brought to the lower resolution of the E(J — K s ) 
map for an unbiased comparison (Section l4.ip . Establishing the precise resolution of color excess maps derived from 
2MASS data has not been a high priority and so we have done our own explorations. A fruitful approach was to use 
the power spectrum of the map. A traditional map (image) produced by an instrument with a 'beam' (or point spread 
function) is typically oversampled relative to the beam size and so the information in neighboring pixels is correlated. 
Therefore, there is a predictable roll-over of the power spectrum at high spatial frequencies corresponding to the size 
of the beam. Although the described Gaussian-weighted and truncated sampling used in 'AvMAP' is not exactly 
convolution, nonetheless we found an empirical roll-over of the power spectrum that can be described quantitatively 
by a beam of FWHM 3(5, which seems reasonable. 

A. 3. Sampling the Herschel T1200 map 

Ideally, the two independent measures of column density, E(J — K s ) and T1200, would sample the sky in the same 
way. Rather than trying to bring the higher resolution T1200 map to the lower effective resolution of the E{ J — K s ) 
map by smoothing^ a more direct approach is to sample the T1200 map at the positions of 2MASS stars, forming the 



15 In preliminary work, the Herschel maps of T1200 were 
smoothed and then regridded. For smoothing, we tried both a 
simple boxcar average of 9 by 9 pixels (11'/ 5 pixels) and convolu- 
tion with a 2(5 Gaussian, finding very similar results. Compared 



to these degraded T1200 maps, our 'AvMAP'-produ ced maps of 
E(J — K s ) (see also Figure 1 of Lombardi et al. 2011) are not ob- 
viously more blurry, i.e., missing all the small-scale high column 
density features, but this assessment is qualitative. 
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weighted average in the same way as in 'AvMAP.' We found that simply extracting cataloged 2MASS stars according 
to the criterion of photometric S/N > 7 produced a weighted number of stars close to 'AvMAP.' While this is not 
precisely an identical sampling, nor is the original T1200 known to pencil-beam precision. Thus we obtain a map of 
T1200 on the same pixels as for E{ J — K s ) and with comparable sampling of the column density and resolution pixel by 
pixel. Quantitatively, there is a roll-over in the power spectrum of the sampled r 12 oo map for each of the two Herschel 
fields that confirms a typical resolution 3'5. 

A. 4. Calibration of the Column Density 
iMartin et al.l (|2012D made a direct calibration of A/h in terms of E{ J — K s ): 

N R = (11.5 ± 0.5) x 10 21 E(J - K s ) cm" 2 . (Al) 

We adopted this empirical relation even though it is calibrated up to only Nh ~5x 10 21 cm -2 . The column density 
in this field, gauged in the near-infrared, ranges to almost an order of magnitude larger, but such column densities are 
not accessible in the calibration because that depends on measurements of H and H 2 column densities made in the 
ultraviolet. 

Instead of A/h, column density from these 2MASS-dcrived maps is often expressed in terms of visual extinction, Ay. 
We prefer not to do so, because we suspect that Ay/E(J — K s ) will vary in the most dense molecular regions, more so 
than Nn/E(J — K s ). Nevertheless, since this is such standard practice, for comparison with the literature we provide 
Ay as well, using the standard diffuse ISM conversion factor; this Ay is just not to be taken literally. 



A. 5. Effect of Finite Bandpasses on the Color Excess 

Here we evaluate the effect of the finite width of 2MASS photometric filter response (the RSRF: ICohen et al.l l2003) 
on the observed extinction. The attenuated signal for a background star measured through a broad photometric filter 
depends upon various parameters such as the intrinsic SED of the source, the width of the BSBF, the shape of the 
extinction curve, and the column of material along the line of sight. Mathematically, the signal can be expressed as 

iW(r ) = F Vo £ (j^-J * e- r °W»°) a RSRF(v)dv, (A2) 

where the optical depth is tq at the central frequency v§ of the bandpass, a is the exponent of the power-law approxi- 
mation to the infrared extinction, 7 is the spectral index of the background star across the bandpass, and the limits of 
the integration span the bandpass. The apparent extinction is therefore Agitcr = — 2.5 log(i ? git er (To)/i'fiiter(Q))- This 
value can be compared to the 'monochromatic' extinction 1.086 t . The J band is most affected because of its higher 
extinction. 

Specifically, we investigated the effect on the color excess E(J — K s ) which is the difference between extinction Aj 
and Ak s - Figure [TOl shows the dependence of the ratio of the apparent color excess E(J — K s ) and the m onochromatic 
color excess on iVjj, where iVjj is obtained from the monochromatic color excess using Equation (|A1[) . Even for the 
lines of sight in Orion A with the highest column density, A/h ~ 50 x 10 21 cm~ 2 , the apparent color excess (and the 
column density that would be deduced from it) is lower by only 5%. For this calculation, we used a = 1.8 and 7 = 1, 
but the results for these column densities are not sensitive to these choices. 

B. PREDICTED 100 fiu EMISSION 

From the SED fits pixel by pixel to the four longest wavelength Herschel bands (Section |3|), we have predicted the 
100 /im brightness and then convolved and re gridded this predicted map for comp arison with the 100 firm IRIS image 
(Improved Reprocessing of the IRAS Survey; iMiville-Deschenes fc Lagachi 2005)13 We compared the maps in two 
ways. First, we computed the power spectrum of each. The shapes of the respective power spectra were very similar. 
Their relative amplitude indicated that the predicted map was 0.93 times fainter on average, very close considering 
that we made no color corrections for the finite IRAS bandpass. 

Second, we made a scatter plot for a pixel by pixel comparison of the brightness in each map. Figure 1111 shows that 
the correlation is very good. For a formal linear fit we used the IDL routine SIXLIN (Iso be et al.lll990T) . adopting the 
bisector. The slope 0.93±0.01 is the same as obtained from the power spectra and the intercept — 1.5±0.5 MJy sr _1 is 
close to zero as expected. The deviation of data from the correlation line at low surface brightness (< 25 MJy sr -1 ) is 
probably related to the uncertainty associated with the zero-point offsets applied to the Herschel maps (Appendix [Cj . 
This comparison validates the SED parameters obtained using only Herschel data and is consistent with the expectation 
that most of the emission at 100 /jm is due to thermal emission from big grains. 

We carried out a similar exercise for 70 /jm and 60 fim emission and found that the brightness was significantly 
underpredicted, as anticipated because of the contribution of non-equilibrium emission by smaller grains. The PACS 
datum plotted in Figure [T] illustrates this point. 



16 Alternatively, we convolved and regridded the SPIRE and . , . , ,. , c 4 . ,. ,„„ 

„.„„ . , J ' TT -, T „ , j ■ °., nj „„„ . , , pixel at this coarser resolution, beiorc predicting the 100 urn emis- 

PACS images to the IRIS resolution and then tit SEDs pixel by . . .. , r j 

a sion; no systematic change was iound. 
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Figure 10. Ratio of synthetic color-excess E(J — K B ) obtained using 2MASS filters to the corresponding monochromatic color excess, as 
a function of Nn. The vertical line shows the upper limit to TVjj in the present analysis of lines of sight in the Orion A molecular cloud, 
and so the effect on this study is small. 

C. ASSESSMENT OF UNCERTAINTIES IN THE SED PARAMETERS 



Errors on the parameters from the SED fit were estimated by Monte Carlo simulation (Chapin et al. 2008). For this 
assessment, we needed to estimate various sources of error in the intensities of the Herschel images. If these errors are 
realistic, then a typical SED fit would have a reduced \ 2 close to unity and a map of \ 2 would be featureless. 

We derived the minimum absolute error for the intensities in each map from the power spectrum analysis as described 
in Section [5] It was clear that a constant error across the map of this order was too small in areas of high brightness; 
for the latter, a constant percentage error of order 5% was required to produce the expected 'flat' behavior in the % 2 
map. 

While investigating the details of the fits of individual SEDs for the Orion A-Cl field we found that the 350 and 
500 /im data points were consistently above and below the best-fit SEDs for surface brightnesses less than 25 and 10 
MJy sr _1 , respectively, and the systematic deviation increased with decreasing surface brightness. This motivated 
us to refine the offsets in the 350 and 500 /jm maps by correlating with the 250 /zm data. For the Orion A-Cl field 
the correlation line between 250 and 160 /im passed through the origin as expected. The correlation lines of 350 and 
500 /im relative to 250 /im had intercepts of 3 and —1.5 MJy sr _1 , respectively. These values are of the magnitude and 
sign that would produce the observed systematic deviations from the best-fit SEDs for low surface brightness pixels. 
The two long-wavelength maps were therefore corrected by an additive offset, making the intercepts formally zero. A 
similar exercise on the Orion A-Sl field also led us to correct the 350, 500, and as well as 160 /jm images by 2.5, 1, 
and —12.0 MJy sr" 1 , respectively. These refined offsets are not perfect and we estimate an error approximately 10% 
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Figure 11. Correlation between the predicted 100 (im intensity obtained by extrapolating the fitted SED and the intensity directly from 
the 100 /an IRIS image. The slope 0.93 ± 0.01 of the correlation line is close to that anticipated from comparison of the power spectra. 
Pixels representing two-dimensional histogram is same as described in Figure [3] Individual data points are plotted where the density of 
points is less than or equal to three. 

of the originally applied zero-point offsets from Planck and IRAS. 

Finally, we estimated the total error in the intensity to be fit by adding in quadrature the minimum absolute error, 
10% of the offset value, and 5% of the intensity. These estimates ought to be conservative and do produce a flat % 2 
image near two (as expected for four intensities, two parameters, and thus two degrees of freedom). Note that the 
actual values of the derived parameters are not sensitive to the precise details of the estimate of total noise. 

In the Monte Carlo simulation of the uncertainties, 500 realizations of mock data were generated by adding inde- 
pendent Gaussian noise consistent with the above values to the actual intensities. For each realization, the SED was 
fit and the corresponding parameters were recorded. An uncertainty on each parameter was obtained by fitting a 
Gaussian to the histogram of the generated distribution. By keeping a record of each fit, we also kept track of the 
correlations of the uncertainties and so can produce the elliptical \-a confidence intervals in, for example, the T - T1200 
plane (see Figure [5]). 
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